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helps  plan  over  $  1  trillion.  Several  approaches  have  been  proposed  and  implemented  to 
solve  the  CIPA  core,  a  mixed-integer  linear  program  (MILP).  Unfortunately,  some  of 
these  MILPs  cannot  be  solved  in  a  reasonable  amount  of  time  using  general-purpose 
commercially  available  optimization  software.  This  thesis  presents  a  new  MILP-based 
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EXECUTIVE  SUMMARY 


Capital  Investment  Planning  Aid  (CIPA)  is  an  optimization  based  decision 
support  system  for  the  U.S.  Navy.  CIPA  prescribes  annual  force  structure  procurement 
and  retirement  plans  based  on  industrial  and  budget  constraints,  as  well  as  mission 
inventory  and  force  mix  requirements  over  a  30-year  horizon. 

CIPA’s  core  is  a  mixed-integer  linear  problem  (MILP).  CIPA  currently  uses  two 
methods  for  solving  the  MILP:  A  heuristic  solver  (HS)  and  an  exact  solver  (ES)  based  on 
Branch-and -Bound  (B&B).  The  HS  is  a  customized  local-search  heuristic.  The  ES  is 
general-purpose  commercial  solver.  Unfortunately,  some  CIPA  instances  cannot  be 
solved  optimally  by  the  HS,  or  feasibly  with  the  ES,  in  reasonable  time.  This  thesis 
presents  a  new  MILP-based  heuristic  technique,  fix-and-relax  (F&R),  which  yields  good 
quality  solutions  and  reduces  the  ES  computational  solution  time  for  our  set  of  realistic 
test  cases. 

This  thesis  also  analyzes  the  “set  aside  budget”  and  “set  aside  labor”  concepts 
implemented  in  CIPA  to  diminish  the  impact  of  end-effects.  The  result  is  a  more  realistic 
problem,  where  procurement  levels  are  better  aligned  with  mission  requirements.  We 
show  the  use  of  these  concepts  mitigates  over-expenditures  in  our  test  cases. 

As  a  result  of  this  research,  we  recommend  using  the  existing  heuristic  solver 
first,  along  with  any  exact  lower  bound.  If  this  does  not  provide  a  solution  within  the 
desired  tolerance,  then  F&R  should  be  used  as  a  second  solver.  If  F&R  cannot  find  a 
solution  within  tolerance,  then  B&B  can  be  used. 
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I.  INTRODUCTION 


“Where  are  the  carriers?”  The  nation's  leaders  ask  that  question  whenever  a  crisis 
involving  vital  U.S.  interests  develops.  Not  only  carriers,  but  also  other  naval  assets  such 
as  destroyers,  frigates,  mine  hunters,  landing  ships,  and  submarines  are  necessary  for 
homeland  security.  Navy  forces  have  a  legitimate  role  in  the  execution  of  national 
security  strategy,  but  like  other  military  systems  can  be  a  target  of  public  and 
congressional  criticism  because  of  their  enormous  cost  [Isenberg,  2002], 

The  construction,  deployment,  and  use  of  naval  assets  have  political  as  well  as 
financial  costs.  Thus,  it  is  important  to  determine  whether  there  are  less  costly  ways  of 
planning  the  procurement  and  retirement  of  naval  assets  to  carry  out  the  missions 
assigned  to  them.  Capital  Investment  Planning  Aid  (CIPA)  has  been  developed  to  help 
navy  analysts  plan  the  retirement  and  procurement  schedules  of  Navy  assets  over  a  30- 
year  planning  horizon.  The  CIPA  core  is  a  mixed-integer  linear  program  (MILP). 
Unfortunately,  some  of  these  MILPs  cannot  be  solved  in  a  reasonable  time  [Salmeron  et 
al.,  2002],  This  thesis  offers  a  new  approach  for  solving  the  existing  CIPA  MILP  that 
achieves  accurate  results  in  most  of  the  existing  test  cases.  The  thesis  also  covers  a 
potential  extension  of  the  model  to  deal  with  end-effects. 

A.  U.S.  NAVY  BUDGET  PLANNING 

Aircraft  carriers  are  a  central  part  of  U.S.  defense  strategy.  The  current  estimate  to 
develop  and  build  the  first  CVNX,  the  next  generation  U.S.  Navy  aircraft  carrier,  is  over 
$10  billion  [Scarborough,  2002],  How  is  this  defense  budget  balanced  amongst  other 
investments  made  by  the  Department  of  Defense  (DoD)? 

DoD  uses  a  Planning,  Programming,  and  Budgeting  System  (PPBS)  to  map  the 
best  course  of  action  to  accomplish  its  missions.  PPBS  is  a  formal,  systematic  structure 
for  making  decisions  on  policy,  strategy,  and  the  development  of  forces  and  capabilities 
to  accomplish  anticipated  missions.  PPBS  is  a  cyclic  process  containing  three  distinct,  but 
interrelated  phases:  planning,  which  produces  defense  planning  guidance;  programming, 
which  produces  approved  program  objectives  memorandum  for  the  military  departments 
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and  defense  agencies;  and  budgeting,  which  produces  the  DoD  portion  of  the  President's 
national  budget  [The  Chief  of  Naval  Operations  (CNO)  N6,  2002]. 

Integrated  Warfare  Architecture  Assessment  and  Planning  Process  (IWARS)  is 
part  of  the  U.S.  Navy  planning  process.  IWARS  comprises  five  warfare  areas  (Power 
Projection,  Sea  Dominance,  Air  Dominance,  Infonnation  Superiority/Sensors, 
Deterrence)  and  seven  support  areas  (Sustainment,  Infrastructure,  Manpower  & 
Personnel,  Readiness,  Training  &  Education,  Technology,  Force  Structure),  which  reflect 
the  complexity  of  naval  warfare  requirements  and  the  need  to  integrate  them  fully  with 
careful  allocation  of  scarce  resources.  Each  of  the  12  IWARS  is  assessed  in  an  attempt  to 
answer  the  question  of  “how  much  is  enough?”,  both  in  terms  of  quality  and  quantity, 
today  and  in  the  future  [CNO  N6,  2002], 

B.  FORCE  STRUCTURE  ANALYSIS 

CIPA  recommends  the  best  yearly  force  structure  procurement  plan  that  satisfies 
industrial  and  budget  constraints  as  well  as  mission  inventory  and  force  requirements. 

CIPA  has  two  methods  to  solve  its  core  MILP,  a  heuristic  solver  (HS)  and  an 
exact  solver  (ES)  [Sahneron  et  ah,  2002]. 

The  HS  is  a  customized  local  search  heuristic  that  typically  returns  a  plan 
satisfying  the  specific  requirements  in  a  matter  of  seconds.  The  solution  accuracy, 
however,  is  case-dependent.  The  HS  also  yields  a  valid  lower  bound,  which  can  be  used 
as  an  objective  assessment  of  the  worst-case  quality  of  the  solution  returned. 

The  ES  attempts  to  solve  the  MILP  exactly.  CIPA  uses  The  General  Algebraic 
Modeling  System  (GAMS)  [Brooke  et  ah,  1998],  a  commercial  algebraic  modeling 
language  to  generate  the  MILP,  and  solves  it  using  a  contemporary  commercial  solver 
(e.g.,  OSL  [GAMS-OSL,  2002],  CPLEX  [GAMS-CPLEX,  2002]).  Unfortunately,  some 
CIPA  instances  cannot  be  solved  optimally  in  a  reasonable  time  using  this  general- 
purpose  optimization  software.  This  thesis  presents  a  new  MILP-based  heuristic 
technique,  fix  and  relax  (F&R),  which  yields  faster  answers  than  ES  without 
compromising  the  quality  of  the  solution  obtained  for  all  test  cases  considered. 
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C.  STATE  OF  THE  ART 

Three  theses  and  one  report  have  been  published  to  date  about  CIPA. 

1.  Planning  Capital  Investments  in  Navy  Forces 

Field  [1999]  presents  the  first  integer-linear  program  of  CIPA.  Field  tests  CIPA 
using  a  2  5 -year  planning  horizon  with  eight  mission  areas,  19  ship  classes,  five  aircraft 
types,  five  production  facilities,  and  three  categories  of  money. 

2.  Optimizing  Procurement  Planning  of  Navy  Ships  and  Aircraft 

Baran  [2000]  introduces  Generalizing  Procurement  Planning  for  Naval  Ships  and 

Aircraft  (GENSA),  which  extends  the  previous  version  of  CIPA.  GENSA  is  tested  with  a 
30-year  planning  horizon  with  29  mission  areas,  45  ship  classes,  39  aircraft  types,  13 
production  facilities,  and  four  categories  of  money. 

3.  Optimized  Procurement  and  Retirement  Planning  of  Navy  Ships  and 
Aircraft 

Garcia  [2001]  focuses  on  improving  the  underlying  optimization  modeling  for 
aircraft  procurement  and  retirement  scheduling  in  the  Capital  Investment  Planning  Aid 
with  Air  Planning  Update  (CIPA  APU).  CIPA  APU  explicitly  incorporates  the  increase 
in  operations  and  maintenance  (O&M)  costs  of  an  aircraft  by  age,  and  deals  with 
retirement  issues  by  aircraft  type  and  age  rather  than  simply  by  aircraft  type. 

4.  Capital  Investment  Planning  Aid  (CIPA)  -an  Optimization-Based 
Decision-Support  Tool  to  Plan  Procurement  and  Retirement  of  Naval 
Platforms 

Sahneron  et  al.  [2002]  describe  the  planning  environment  into  which  CIPA  has 
been  introduced,  showing  how  CIPA  works,  and  how  CIPA  is  used.  The  report  presents 
an  overview  of  CIPA.  It  describes  the  planning  environment,  and  presents  the  latest 
version  of  the  underlying  MILP  at  the  heart  of  CIPA,  discusses  exact  and  heuristic 
techniques  used  to  solve  CIPA,  along  with  their  computational  performance,  and 
provides  an  overview  of  the  graphical  user  interface.  Since  this  report  is  the  latest 
research  about  CIPA,  it  contains  the  latest  version  of  the  MILP  used  in  this  thesis. 
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II.  CIPA  MIXED-INTEGER  LINEAR  PROGRAM 


A.  MATHEMATICAL  MODEL  OVERVIEW 

CIPA  minimizes  penalties  associated  with  violating  budget  constraints, 
production  constraints,  and  inventory  requirements.  CIPA  gives  a  recommended  plan  that 
includes  budget,  purchase  dates,  quantities  and  cost,  and  production  facility  work-force 
levels.  CIPA  also  isolates  force  level  deficiencies  inflicted  by  budget  restrictions  on 
procurements,  production  that  cannot  keep  pace  with  procurement  requirements,  or  the 
lack  of  any  existing  replacement  for  retired  platforms.  CIPA  maintains  yearly  time 
resolution  for  25  or  30  years.  Since  it  can  take  up  to  nine  years  to  build  platforms  such  as 
destroyers,  frigates,  and  submarines,  CIPA’s  prescriptions  for  the  last  few  years  of  the 
planning  horizon  may  suffer  from  end  effects.  The  solution  for  the  last  years  of  the 
horizon  may  not  be  accurate  because  information  for  years  beyond  the  horizon  has  not 
been  specified. 

The  MILP  represents  a  number  of  features  divided  into  six  categories  [Salmeron 
et  ah,  2002]. 

1.  Mission: 

•  Ship-mission  and  air-mission  requirements 

2.  Inventory: 

•  Initial  inventory  of  ships  and  aircraft 

•  Ongoing  (resident)  production  of  ships  and  aircraft 

•  Minimum  and  maximum  annual  production  of  ships  and  aircraft 

•  Maximum  total  production  of  ships  and  aircraft 

•  Maximum  annual  inventory  of  ships  and  aircraft 

•  Minimum  and  maximum  annual  ship  and  aircraft  retirement 

3.  Cost: 

•  Ship  and  aircraft  cost  profile 

•  Economy-of-scale  for  ship  and  aircraft  procurement 

•  Operation  and  maintenance  costs  for  each  ship  and  aircraft 
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4.  Budget: 

•  Minimum  and  maximum  annual  budget  available 

•  Minimum  and  maximum  cumulative  budget  available 

•  Set  aside  budget  (for  ships  and  aircraft) 

5.  Industry: 

•  Work- force  profile  for  ship  production 

•  Minimum  and  maximum  annual  work-force  levels  for  ship 
industry 

•  Set  aside  labor  for  ships  and  its  relationship  with  set  aside  budget 

6.  Penalty: 

•  Tradeoff  among  budget  shortfall  (or  surplus),  industry  work-force 
shortfall  (or  surplus)  and  mission  shortfall 

Mission  requirements  (category  1)  drive  platform  procurement.  Category  2 
features  account  for  yearly  platfonn  inventory  levels  as  well  as  shipyard  capacity, 
minimum  retirement  levels  and  the  age  of  existing  platfonns.  Category  3  considers  CIPA 
cost-related  features.  Procurement  costs  are  typically  incurred  and  spread  out  over  a 
number  of  years  before  a  platform  is  delivered.  The  cost  of  purchasing  platforms  exhibits 
economies  of  scale.  Category  4  specifies  annual  and  cumulative  expenditures  that  should 
not  exceed  or  fall  below  their  respective  specified  limits.  Category  5  refers  to  work-force 
requirements  for  ship  production  that  are  spread  out  over  the  production  period  of  a  ship. 
Ideally,  workforce  levels  should  stay  within  specified  limits  to  prevent  the  loss  of 
industrial  capability  and  to  avoid  overtime  costs.  The  last  category  refers  to  CIPA  penalty 
charges  for  each  individual  violation  of  budget,  industry,  or  mission-required  levels.  The 
penalties  express  the  tradeoff  among  the  different  shortfalls  and  surpluses  in  order  to 
prioritize  the  satisfaction  of  those  conditions  deemed  more  critical  by  the  planner. 

As  main  decision  variables,  the  number  of  platforms  procured  and  retired  every 
year  is  considered.  Additional  variables  are  added  to  specify  the  piece-wise  linear 
approximation  of  non-convex  costs  associated  with  economies-of-scale.  “Elastic” 
variables  are  also  incorporated  to  account  for  budget,  industry,  and  mission  requirement 
violations.  The  objective  function  minimizes  the  sum  of  these  violations.  See  Field 
[1999]  for  a  discussion  of  how  to  select  penalty  values. 
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All  these  features  are  mathematically  represented  through  the  following  linear 
program: 

CIPA:  min  F 

s.t.  (1)  to  (50) 

where  the  objective  function,  F,  and  the  constraints  (1)  and  (50),  are  described  in  detail  in 
the  following  section. 

B.  CIPA  MODEL 

This  section  presents  the  mathematical  fonnulation  of  the  CIPA  model  presented 
in  Sahneron  et  al.  [2002],  page  17,  and  incorporates  the  proposed  changes  for  end-effects 
described  on  page  1 19  of  the  same  document.  The  formulation  of  the  model  is  included  in 
this  document  to  be  comprehensive. 

1.  Sets  and  Indices 

•  Time 

Y,  set  of  years  of  the  planning  horizon;  y,y'e  Y  .  For  convenience,  it 

is  assumed  that  Y  =  {1,2,3..., |  Y  |} 

•  Platform 

A,  set  of  aircraft  types;  a  e  A 

S,  set  of  ship  classes;  se  S 

•  Mission 

M  A ,  set  of  air  missions;  m  e  M  A 

Ms ,  set  of  ship  missions;  m  e  Ms 

Am  cz  A ,  subset  of  aircraft  types  that  contribute  to  mission  m  e  MA 

Sm  c  S  ,  subset  of  ship  classes  that  contribute  to  mission  m  e  Ms 

•  Production 

Ia ,  set  of  cost  increments  for  aircraft  a  e  A;  ie  Ia 

P,  set  of  production  facilities;  pe  P 

Ps  ci5  ,  subset  of  facilities  that  produce  ship  class  s  e  S 


7 


Q  ,  set  of  quantities  available  for  ship  se  S  procurement  at  facility 

p  e  Ps  in  year  y  e  Y .  This  set  is  defined  in  tenns  of  the  sproc 

and  sproc  parameters  (see  below)  as  follows: 

q  e  Qspy  =  {sproc  ^  , sproc _  +!,■■■, sproc spy } 

•  Others 

Z  + ,  set  of  non-negative  integers,  Z+  =  {0,1,2,...} 

2.  Parameters  (and  Units) 

•  Conventions 

The  word  “procurement”  or  “to  procure”  refers  to  “delivery”  or  “to  deliver”, 
respectively,  unless  explicitly  stated  otherwise.  Therefore,  “procure”  is  referred  to  as  the 
action  that  takes  place  at  the  moment  (year)  that  the  platform  is  delivered  and  available 
for  use  from  that  year  onwards,  regardless  when  the  real  “procurement”  arrangements 
were  made. 

The  words  “time  period”  and  “year”  are  used  interchangeably. 

The  words  “facility”  and  “plant”  are  used  interchangeably. 

•  Objective-related  parameters:  Penalties 

ampenm,  penalty  for  shortage  in  completing  air  mission  m  e  M  A  ( $  per 
aircraft) 

smpenm,  penalty  for  shortage  in  completing  ship  mission  me  Ms  ($  per 
ship) 

bpen+ ,  penalty  for  budget  excess  ($  per  $) 

bpen~ ,  penalty  for  budget  shortage  ($  per  $) 

cbperiy  ,  penalty  for  cumulative  expenses  excess  ($  per  $) 

cbpen~ ,  penalty  for  cumulative  expenses  shortage  ($  per  $) 

lpen+p ,  penalty  for  labor  excess  at  plant  p  e  P  ($  per  worker) 

lpen~ ,  penalty  for  labor  shortage  at  plant  p  e  P  ($  per  worker) 
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•  Constraint-related  parameters:  used  for  indices  dependencies 

SBb  ,  number  of  years  before  (starting  at  0)  the  procurement  of  ship  class 

se  S  from  plant  p  e  Ps  requires  budget  (i.e.,  in  0,1,...  SBbsp  -1 
years  before) 

SCbsp,  number  of  years  before  (starting  at  0)  the  procurement  of  ship  class 

se  S  from  plant  pe  Ps  requires  labor  (i.e.,  in  0,1,...  SCb  -1 
years  before) 

SBasp,  number  of  years  after  (starting  at  1)  the  procurement  of  ship  class 

se  S  from  plant  pe  Ps  requires  budget  (i.e.,  in  0,1,...  SBa  years 
before) 

SCasp,  number  of  years  after  (starting  at  1)  the  procurement  of  ship  class 

se  S  from  plant  pe  Ps  requires  labor  (i.e.,  in  0,1,...  SCa  years 
before) 

ABba ,  number  of  years  before  the  procurement  of  aircraft  type  a  e  A  in 

which  the  aircraft  is  paid  (at  once) 

•  Constraint-related  parameters:  Ships 

sin  v°,  initial  inventory  of  class  s  e  S  ships  (number  of  ships) 

csprocsy,  committed  procurement  of  class  se  S  ships  in  year  yeY  due  to 

production  in  progress  (number  of  ships) 

sin vs,  maximum  number  of  class  se  S  ships  in  inventory  (number  of 

ships) 

stotsp,  maximum  number  of  class  se  S  ships  to  procure  from  plant 

p  e  P  (number  of  ships) 

sproc  ,  minimum  number  of  class  se  S  ships  to  procure  from  plant 

p  e  P  in  time  period  ye  Y  (number  of  ships) 

Note:  sproc  =  0 ,\/seS,pePs ;  Vy  <  max {SBbsp ,  SCbsp  }  - 1  and 

sproc  =0,  \/se  S,pe  Ps ;  Vy  > \  Y  \  +1  -  max  {SBasp ,  SCasp  }  is 
required 

sproc  ,  maximum  number  of  class  seS  ships  to  procure  from  plant 

pe  Ps  in  time  period  ye  Y  (number  of  ships) 
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and 


Note:  sprocspy  =  0,  Vs  e  S,  p  e  Ps ;  Vy  <  max  {SBbsp ,  SCbsp  }  - 1 

sproc  =0,  Vs  e  S,  p  e  Ps ;  Vy  > |  Y  |  +1  -  max  {SBasp ,  SCa  }  is  required. 

•  Constraint-related  parameters:  Aircraft 

ainv°a,  initial  inventory  of  type  ae  A  aircraft  (number  of  aircraft) 

caproc  ,  committed  procurement  of  type  a  e  A  aircraft  in  year  ye  Y  due  to 
production  in  progress  (number  of  aircraft) 
ainva,  maximum  number  of  type  ae  A  aircraft  in  inventory  (number  of 

aircraft) 

atota,  maximum  number  of  type  ae  A  aircraft  to  procure  (number  of 

aircraft) 

aproc  ,  minimum  number  of  type  ae  A  aircraft  to  procure  in  time  period 
ye  Y  (number  of  ships) 

aproc  av,  maximum  number  of  type  ae  A  aircraft  to  procure  in  time  period 

ye  Y  (number  of  ships) 

inc  j,  increment  i  e  I  a  lower  bound  for  the  number  of  type  ae  A  aircraft 

to  be  procured  in  year  ye  Y  (number  of  aircraft) 
incayi,  increment  i  e  I  a  upper  bound  for  the  number  of  type  ae  A  aircraft 

to  be  procured  in  year  ye  Y  (number  of  aircraft) 
squad  a ,  squadron  size  for  aircraft  ae  A  procurement  (number  of  aircraft) 

•  Constraint-related  parameters:  Retirements 

csret ,  minimum  cumulative  number  of  class  se  S  ships  to  retire  by  the 

end  of  time  period  ye  Y  (number  of  ships) 
csret sy,  maximum  cumulative  number  of  class  se  S  ships  to  retire  by  the 

end  of  time  period  ye  Y  (number  of  ships) 

sret  ,  minimum  number  of  class  se  S  ships  to  retire  by  the  end  of  time 

period  ye  Y  (number  of  ships) 

sret  sy ,  maximum  number  of  class  s  e  S  ships  to  retire  by  the  end  of  time 

period  vg  Y  (number  of  ships) 

caret  ,  minimum  cumulative  number  of  type  ae  A  aircraft  to  retire  by  the 
end  of  time  period  vg  Y  (number  of  aircraft) 
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caret ay,  maximum  cumulative  number  of  type  ae  A  aircraft  to  retire  by 

the  end  of  time  period  ye  Y  (number  of  aircraft) 

qret  ,  minimum  number  of  type  ae  A  aircraft  to  retire  by  the  end  of 
time  period  ye  Y  (number  of  aircraft) 

aretsy,  maximum  number  of  type  ae  A  aircraft  to  retire  by  the  end  of 
time  period  ye  Y  (number  of  aircraft) 

Constraint-related  parameters:  Missions 


seff  sm , 

effectiveness  for  ship  5  e  S m  performing  mission 
(number  of  missions  per  ship) 

m  e  M  s 

aeffam  > 

effectiveness  for  aircraft  a  e  A tn  performing 

in  e  M  A  (number  of  missions  per  aircraft) 

mission 

smreq  , 

my 

overall  effectiveness  required  for  ship  mission  meMs 
period  ye  Y  (number  of  missions) 

in  time 

amreq  , 

my 

overall  effectiveness  required  for  air  mission  me  MA 
period  ye  Y  (number  of  missions) 

in  time 

•  Constraint-related  parameters:  Budget 
oscn  ,  fixed  SCN  cost  in  year  ye  Y  ($) 

ocscny,  fixed  SCN  cost  in  year  ye  Y  for  ships  not  considered  ($) 

frac,  historical  fraction  of  total  SCN  cost  for  ship  outfitting 

oapny,  fixed  APN  cost  in  year  ye  Y  ($) 

ocapn  ,  fixed  APN  cost  in  year  ye  Y  for  aircraft  not  considered  ($) 

apn5,  historical  fraction  of  total  APN  categories  1  through  4  required  for 

categories  5  through  7 

oomy,  fixed  O&M  cost  in  year  ye  Y  for  maintenance  not  considered  ($) 

scostb  lt  SCN  cost  incurred  /  years  before  q  class-5  ships  are  procured  from 
plant  p,  for  s  e  S ,  pePs,  qe\^  Qspy  ,  /  =  {0,1  — ,SBbsp  - 1}  ($) 

jpey 

s costa  SCN  cost  incurred  /  years  after  q  class-5  ships  are  procured  from 

plant  p,  for  s  e  S,  pePs,  g  e  |J  Qspy  ,  l  =  {l,---,SBasp}  ($) 

jpey 

aacostayi,  increment  ie  I a  procurement  cost  for  type  ae  A  aircraft  in  year 
ye  Y  ($  per  aircraft) 
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abcostayi, 

omship  , 
omairav, 
csbudgety , 

toav, 

toay, 

ctoay, 

ctoay, 

ssaby, 

asaby , 


increment  i e  I a  fixed  procurement  cost  (intercept)  for  type  ae  A 
aircraft  in  year  ye  Y  ($) 

O&M  cost  for  class  se  S  ship  in  year  ye  Y  ($  per  ship) 

O&M  cost  for  type  ae  A  aircraft  in  year  ye  Y  ($  per  ship) 
committed  budget  in  year  vg  Y  due  to  ship  production  in  progress 
($) 

TO  A  budget  lower  limit  for  year  ye  Y  ($) 

TO  A  budget  upper  limit  for  year  ye  Y  ($) 

TO  A  cumulative  budget  lower  limit  for  year  ye  Y  ($) 

TO  A  cumulative  budget  upper  limit  for  year  ye  Y  ($) 
maximum  set  aside  ship  budget  for  year  vg  Y  ($) 
maximum  set  aside  aircraft  budget  for  year  vg  Y  ($) 


Constraint-related  parameters:  Labor 


claborpy, 

sworkbspqn, 


sworkaspqn, 


pcap 


py’ 


PCCiPpy’ 


sal  py, 


lcratep, 


committed  labor  in  year  vg  Y  at  plant  pe  P  due  to  production  in 
progress  (number  of  workers) 

required  labor  n  years  before  q  class-5  ships  are  procured  from 
plant  p,  for  seS,  pePs,  qe{]Qspy,  n  =  {Q\,---,SCbsp -1} 

yeY 

(number  of  workers) 

required  labor  n  years  after  q  class-5  ships  are  procured  from  plant 
p,  for  seS,  pePs,  q^\jQspy,  n  =  {1, •  - ■  ,SCasp)  (number  of 

yeY 

workers) 

minimum  production  capacity  at  plant  pe  P  in  time  period  vg  Y 
(number  of  workers) 

maximum  production  capacity  at  plant  pe  P  in  time  period  vg  Y 
(number  of  workers) 

maximum  set  aside  labor  at  plant  pe  P  in  time  period  ye  Y 
(number  of  workers) 

approximate  labor  cost  at  plant  pe  P  for  set  aside  labor  purposes 
($/worker) 


3.  Decision  Variables  (and  Units) 

•  Variables  related  to  objective  function  and  to  elastic  constraints 
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objective  function  value 

air  mission  m  e  M A  shortage  in  year  ye  Y  (number  of  aircraft) 
ship  mission  m  e  Ms  shortage  in  year  ye  7  (number  of  ships) 
budget  excess  in  year  ye  7  ($) 
budget  shortage  in  year  yef  ($) 
cumulative  budget  excess  in  year  yef  ($) 
cumulative  budget  shortage  in  year  yef  ($) 
labor  excess  in  year  ye  7  (number  of  workers) 
labor  shortage  in  year  ye  7  (number  of  workers) 

•  Main  decision  variables 

APROCavi,  number  of  type  ae  A  aircraft  to  procure  at  the  start  of  year  ye  Y 
in  cost  increment  i  e  I a  (number  of  aircraft) 

ARETav,  number  of  type  a  e  A  aircraft  to  retire  by  the  end  of  year 
ye  Y  (number  of  aircraft) 

SPROC  ,  one  if  facility  pe  P  is  to  deliver  qe  Q  class  se  S  ships  at  the 
start  of  year  ye  Y ,  and  zero  otherwise  (0-1  variable) 

SRET '  number  of  class  seS  ships  to  retire  by  the  end  of  year 

ye  Y  (number  of  ships) 

SSABudget  ,  amount  of  budget  set  aside  in  year  ye  Y  for  future  ship 

procurements  ($) 

AS  A  Budget y ,  amount  of  budget  set  aside  in  year  ye  Y  for  future  aircraft 
procurements  ($) 

SALabor  ,  amount  of  labor  set  aside  in  year  yeY  for  future  ship 

procurements  from  plant  pe  P  (number  of  workers) 

•  Control  decision  variables 

APavi,  one  if  aircraft  ae  A  is  procured  at  the  start  of  year  ye  Y  in  cost 

increment  i  e  Ia ,  and  zero  otherwise  (0- 1  variable) 

A1NV  ,  inventory  of  type  ae  A  aircraft  at  the  start  of  year  ye  Y  (number 
of  aircraft) 

AMEff  ,  overall  effectiveness  achieved  for  air  mission  me  MA  in  year 
ye  Y  (number  of  missions) 


F, 

tAM 

my 

SM 
*  my  ' 

tB+ 
ly  ’ 


az 


a  f, 


CB- 
y  ’ 
L+ 

y  ’ 
L- 
y  ’ 
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SINVsy,  inventory  of  class  s  e  S  ships  at  the  start  of  year  ye  Y  (number  of 
ships) 

SMEff  ,  overall  effectiveness  achieved  for  ship  mission  me  M s  in  year 
ye  Y  (number  of  missions) 

SBUDGETy,  amount  of  SCN  money  to  budget  for  year  ye  Y  ($) 

ABUDGETv,  amount  of  APN  money  to  budget  for  year  ye  Y  ($) 

OMBUDGET v,  amount  of  O&M  money  to  budget  for  year  ye  Y  ($) 

BUDGET Y,  total  amount  of  money  to  budget  for  year  ye  Y  ($) 

LABOR  ,  amount  of  labor  required  in  year  ye  Y  at  plant  pe  P  (number  of 
workers) 

4.  Formulation 

mill  F  =  Yj  YJamPenm  °C  +  X  YusmPenm  °C  + 

yeY  meMA  yeY  meMs 

YjbPeK  <  +  YjhPen~y  ay~  +  TjcbPen+y  a7+  +  TjcbPeny  aT~  + 

yeY  yeY  yeY  yeY 

XX  +  XZH  aLpy 

yeY  peP  yeY  peP 

subject  to: 

•  Ship 

YSPR0C.m=  >"  VseS,pePs;VyeY  (1) 

qeQspy 

SINVsy  =  sin  v°  +  I  csprocsy,  +  II  Yl  SPROC.p,;  -  ZSRET.y 

y'eY\y'<y  pePsy'<yqeQspy.  y'eY\y'<y-l 

\/se  S;\/ye  Y  (2) 

X  1ll1  SPR0Cspyq  ^  st°Gp,  V.v  e  S,pe  P  (3) 

y^Y  qeQw 
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Aircraft 


iela 

Mayi  APayi  <  APROCayi  <  incayi  APayi, 
aPr°C„y  ^  YjAPR0Cayi  -  aPrOCay  > 


Vae  A;Vye  Y 


\/ae  A, is  7a;Vye  7 
Vo  e  A;\/ye  Y 


AINVay  =  ainvl  +  £. capro, V  +  £  YAPR0Cy<~  YARETV" 

y'eY\y'<y  y'eY\y'<y  iela  y’eY\y'<y— 1 

Vae  A;Vye  Y 


APROCayi  <  atota,  Vae  A 

yeY  iela 

•  Retirements 


csret_sy  <  ^SRETy,  <  csretsy, 

y'eY\y’<y 


caret  „y  <  ^  ARETay ,  <  caret ay, 

v'eY\y'<y 


\/s  €  S;Vye  Y 
Va<E  A;Vye  Y 


•  Mission  Inventory 

*/,  =  sink,, 

seSm 

SMEffmy  +  *»  >  Jmreg_, 

Am(fm,  =  AINV-r 

as  Am 

AMEff my  +  a ™  >  amreq  , 


Vflie  M';Vve  7 

VfflG  M';Vve  7 
Vme  MA;\/ye  Y 

Vme  MA;\/ye  Y 


(4) 

(5) 

(6) 


(7) 

(8) 


(9) 

(10) 


(11) 

(12) 

(13) 

(14) 
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Budget 


SBUDGETv  =  oscny  +  (1  +  frac)(o  esc  ny  +  cs  budget  y  + 

yy  y  y  jcos^  ,  sproc  ,  + 

/—!  /—!  l—t  L-i  spq,y'-y  spy  q 

seS  pei5,  v'eY\ 

y<y'<v+SBbsp 

\ 

y  y  y  y  ^cos  taspqy_y,  sproc spy,q  , 

seS  pePs  v'eY \  qeO  , 

Vye7  (15) 

ABUDGET 7  =  oapnv  +  (1  +  apn5){ocapn  + 

y  y  (««  cos  APROC a  nABb  i  + 

as  A  is  I  a 

abcostay+ARhJ  APay+ABba  i)), 

Vye7  (16) 

OMBUDGETv  =  oom,  +  y  omshipsv  SINVsv  +  y  omair av  AINVav,  Vye  7  (17) 

BUDGET  =  SBUDGET  +  ABUDGETV  +  OMBUDGETv  + 

y  y  y  y 

SSABudget  y  +  ASABudget  y , 

VvgF  (18) 

toa>;  <  <xj“  +  BUDGET y,  Vve  7  (19) 

BUDGET v  -  aBv+  <  toay,  Vv  e  7  (20) 

croaj.  <  aCT~  +  YjBUDGETv>  Vvg  7  (21) 

y  B  UDGETV,  -  a  f +  <  Doa  v ,  Vv  e  7  (22) 
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Industrial 


LABORpy  =  claborpy  +  SALaborpv  + 


z  z 

5GiS|/?eP9  y'eY\ 

y<v'<y+SCbsp 

Y;SWOrkbsPq,y'-y 

SPROCspy,q  + 

z  z 

seS\pePs  y'eY \ 

y-SCasp<y'<y 

TjsworkasPq,y- 

^  4eQspy' 

-y  SPROC spy,q , 

\/p€  P;Vy  e  Y 

(23) 

Pcap  pv  ^  ttpy  +  LABOR py, 

Vpe  P;Vy  e  Y 

(24) 

LABOR py  -  aLpy  <  pcap py, 

Vpe  P;Vy  e  Y 

(25) 

crate  SALaborpv  =  SSABudget v , 

psP 

Vpe  P;Vy  e  Y 

(26) 

•  Non-negativity  and  bounds 

0  <  SSABudget y  <  ssaby, 

Vye  Y 

(27) 

0  <  ASABudget v  <  asaby, 

Vye  Y 

(28) 

0  <  SALaborpy  <  sal  py, 

Vpe  P;\/ye  Y 

(29) 

0  <  AINVav  <  aim a 

Vae  A;Vye  Y 

(30) 

AMEffmy  >  0, 

Vffle  MA;\/ye  Y 

(31) 

0  <  SINVsv  <  sin  vs , 

\/se  S;\/ye  Y 

(32) 

SMEff  ny  >  0, 

Vflie  Ms;Vy  e  Y 

(33) 

sretsy  <  SRETsv  <  sretsy. 

\/se  S;\/ye  Y 

(34) 

aret  ay  <  ARETav  <  aretay, 

Vne  A;Vye  Y 

(35) 

SBUDGETy  >  0, 

Vye  Y 

(36) 

ABUDGETv  >  0, 

Vye  Y 

(37) 

OMBUDGETv  >  0, 

Vye  Y 

(38) 

BUDGETV  >  0, 

Vye  Y 

(39) 

LABORpy  >  0, 

\fp  g  P;\/y  e  Y 

(40) 
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a  >  0 


(41) 


Fixed  variables 


APROCayi  =  0, 

Vo  G  A,i  G  Ia ;  Vy  G  Y  \  y  <  ABba 

(42) 

SPROCspy0  =  1, 

VseS,  pe  P;  Vy  g  Y\y<  max  {SBbsp ,  SCbJ-l 

(43) 

SPROCspy0  =  1, 

\/s  g  S'^/ig  P  ;  Vy  g  7|  y  >|  7|  +1  -  max  {SBasp ,  SCasp } 

(44) 

•  Binary/Integer  variables 

APROCayie  Z+ , 

Vo g  ^4,/ g  /a;VyG  7  (45) 

ARETayeZ+, 

Vo  g  A;\/ye  Y 

(46) 

APayi*{  0,1}, 

Vog  ^4,/g  /a;VyG  7  (47) 

SPROCspyqe{  0,1}, 

Vs  g  7,/?g  Ps;Vy  g  Y;\/qe  Qspy 

(48) 

SRETyeZ\ 

\/s  g  7;Vy  g  7 

(49) 

An  additional  constraint  requires  that: 

APROCayj  is  a  multiple  of  squad a ,  Vog  A, is  I a ;  Vy g  7  (50) 

Remark:  This  constraint  is  not  explicitly  stated  in  the  formulation.  However, 
notice  that  it  can  be  easily  addressed  by  setting  the  proper  segment  limits.  For  example, 
if  squad a  =  4  then  the  segment  limits  could  be: 

inCgyi  =  0  =  incgyi,incay2  =  4  =  incay 2,  mco;3  =  8  =  incay 3,... 

Notice  that,  unless  squad a  =  1 ,  in  which  case  extra  segments  are  not  needed,  the 

number  of  segments  in  the  model  is  significantly  increased. 

5.  Description  of  the  Formulation 

Specifically,  the  formulation  serves  the  following  purposes: 
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The  objective  function,  F,  comprises  the  sum  of  all  the  penalties  due  to 
Air-Mission  and  Ship-Mission  shortfall,  budget  deficit  and  surplus, 
cumulative  budget  deficit  and  surplus,  and  labor  deficit  and  excess. 

Ship  constraints  (1)  to  (3)  constrain  ship  procurement:  (1)  ensures  that  one 
option  for  ship  procurement  is  executed  yearly  at  each  plant,  (2)  calculates 
the  yearly  ship  inventory,  and  (3)  limits  the  maximum  procurement  from 
each  plant. 

(4)  to  (8)  constrain  aircraft  procurement:  (4)  to  (6)  guarantee  that 
procurements  are  made  within  the  limits  of  one  specific  segment  and 
without  exceeding  the  general  minimum  and  maximum.  (7)  calculates  the 
yearly  aircraft  inventory  and  (8)  limits  the  maximum  total  procurement 
throughout  the  years. 

Cumulative  retirement  goals  are  specified  in  (9)  to  (10). 

(11)  to  (14)  keep  track  of  platform  inventory  to  perform  each  specific 
mission  and  then  calculate  mission  shortfalls,  which  depend  on  the  overall 
effectiveness  achieved  for  each  mission. 

Budget  constraints  (15)  to  (22)  are  as  follows:  (15)  calculates  the  ship- 
budget  per  year,  which  depends  on  the  payment  profile  for  each  specific 
ship  that  has  been  procured,  (16)  is  the  yearly  aircraft  budget,  considering 
the  segment  cost  definition,  (17)  detennines  O&M  costs  based  on  existing 
inventories.  The  total  yearly  budget  is  assessed  in  (18),  which  serves  to 
compute  deficits  and  surpluses  on  a  yearly  and  cumulative  basis  in  (19)  to 
(22).  Notice  that  the  total  budget  computed  in  (18)  includes  the  budget  set 
aside  to  account  for  end-effects. 

Based  on  labor  profiles  for  those  ships  that  have  been  procured,  the  labor 
force  level  required  at  the  different  shipyards  is  estimated  in  equation  (23). 
Then,  the  lack  of  labor  or  excess  is  computed  in  (24)  to  (25). 


•  Accounting  for  end-effects,  (26)  establishes  an  approximate  relationship 
between  set  aside  labor  and  set  aside  budget. 

•  (27)  to  (41)  establish  non-negativity  and  bounds  for  the  decision  variables. 
Among  these  bounds,  specified  maxima  and  minima  for  platform 
inventory  and  retirement  levels  exist,  and  maximum  levels  for  the  set  aside 
budget  and  set  aside  labor. 

•  Some  variables  need  to  be  fixed  in  (42)  to  (44),  since  they  would 
otherwise  involve  actions  beyond  the  horizon  limits. 

•  (45)  to  (49)  specify  those  variables  that  need  to  be  considered  integer  or 
binary.  This  also  implies  the  integrality  of  other  variables  such  as  platform 
inventories  and  mission  inventories. 

•  Finally,  (50)  requires  the  aircraft  procurement  to  be  a  multiple  of  the 
squadron  size.  As  the  remark  indicates,  this  can  be  accomplished  by 
adding  extra  segments  for  those  aircraft  whose  squadron  size  for 
procurement  purposes  is  greater  than  one. 

C.  CONSIDERATIONS  REGARDING  THE  CIPA  MODEL 

The  previous  work  by  Field  [1999],  Baran  [2000],  and  Garcia  [2001]  did  not 

solve: 

CIPA:  min  F 

s.t.  (1)  to  (50) 

but  a  slightly  relaxed  version  we  call  “Simplified  CIPA  (SCIP A)”: 

SCIPA :  min  F 

(1)  to  (44) 

(47)  to  (48) 

s.t.  <  APROCayi  >  0  (45  -  modified) 

ARETay  >  0  (46  -  modified) 

SRETy  >  0  (49  -  modified) 
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SCIPA  relaxes  the  integrality  requirements  for  aircraft  procurement  and 
retirement  and  for  ship  retirement.  Squadron  size  requirements  for  aircraft  procurement 
(50)  are  also  disregarded. 

Other  constraints  such  as  mission  effectiveness  or  set  aside  budget  and  set  aside 
labor  constraints  are  also  not  considered  in  previous  work  by  Field  [1999],  Baran  [2000] 
and  Garcia  [2001]. 

SCIPA  is  used  rather  than  CIPA  because  it  still  provides  helpful  prescriptions  and 
it  should  be  easier  to  solve.  Salmeron  et  al.  [2002]  devised  a  post-processor  that 
heuristically  rounds  a  SCIPA  solution  and  also  satisfies  the  squadron  size  requirements  to 
produce  a  feasible  solution  to  the  original  CIPA  MILP. 

Consequently,  when  referring  to  a  solution  provided  by  the  exact  solver  (ES),  or 
simply,  to  an  “exact  solution,”  the  solution  of  the  following  process  is  actually  being 
referred  to: 

Exact  Solution=  Solve  (exactly)  SCIPA  +  Round  solution  to  meet  (45)-(46)  and  (49)-(50). 
Of  course,  the  “exact  solution”  is,  in  actuality,  a  heuristic  solution  that  by  construction  is 

expected  to  be  relatively  close  to  the  optimal. 

The  thesis  does  not  attempt  to  improve  the  solution  rounding  process.  Instead,  a 
technique  is  devised  to  efficiently  solve  SCIPA. 

D.  SOLVING  THE  SCIPA  MODEL  USING  FIX-AND-RELAX 

In  this  section  we  present  a  general-purpose  technique  intended  to  reduce  the 
computational  burden  of  solving  SCIPA  by  B&B.  The  technique  involves  using  a  number 
of  sub-problems  each  with  fewer  binary  variables  than  the  original  MILP.  The  approach 
follows  Fix-and-Relax  (F&R)  introduced  by  Dillenberger  et  al.  [1994],  See  also  Escudero 
and  Salmeron  [2002]  and  earlier  versions  of  similar  techniques  used  by  Brown  et  al. 
[1987],  among  others. 
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SCIPA  can  be  rewritten  as: 


SCIPA : 


min  f(x,z ) 
ze  {0,1}" 


s.t.l 


x>0 

( x,z)e  XZ 


where  z  is  a  vector  that  comprises  APuyi,\/a  e  A,  i  e  7a;Vv  e  7,  and 
SPROC  ,Vs e  S,  pe  P;\/ye  Y;\/qe  Qspy  .  These  variables  are  required  to  be  binary,  as 

in  the  original  CIPA,  (47)  and  (48).  On  the  other  hand,  x  is  the  vector  of  all  the 
continuous  variables  of  the  SCIPA  model  (APROC,  ARET,  SRET,  a ,  etc.,  with  their 
appropriate  index  domain).  f(x,z )  =  f(ox)  represents  our  linear  objective  function.  The 
set  of  constraints  represented  by  (x,z)e  XZ  is  (1)  to  (44). 


To  generalize  the  exposition  of  the  methodology,  the  components  of  z  are  denoted 
Zj,...,zn,  so  n  is  the  total  number  of  binary  variables  in  the  original  model.  Let 

V  =  {1,2,...,/?}  be  the  set  indices  for  those  variables,  and  let  Vv...,Vk  be  a  direct  partition 

k 

of  the  set  V,  that  is,  Vt  cz  V,\/i  =  =  [Jfj ,  and  Vt  nV,  =  0,V/,/  =  1,...,A:|?  ^  i  . 

i= 1 

The  cardinality  of  each  Vi  is  denoted  |f)  |  =  nt .  Therefore,  n  =  ni ,  and  SCIPA  can  be 


rewritten  as: 


SCIPA :  min  f(x,z ) 

[z,e{0,l},Y/e  V„i  =  l,...,k 


s.t.j 


J 

x>0 


( x,z)e  XZ 


In  the  partition  selected  for  our  problem,  for  a  given  year  y ,  V  is  defined 
comprising  all  the  variables  of  type  APayi,\/ae  A,ie  Ia  and 
SPROC  ,\/s e  S,pe  P;\/qe  Q  ,  i.e.,  all  the  variables  associated  with  year y. 
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The  F&R  framework  solves  the  following  sequence  of  k  mixed-0-1  sub-problems, 
hereafter  called  stages,  and  is  denoted  as  SCIPAr,  for  r  =  In  our  approach, 

k  =  |r| ,  but  in  a  more  general  framework,  k  would  depend  on  how  each  stage  is  defined. 


SCIPAr  is  defined  as  follows: 


SCIPAr 


min  f(x,z) 

( x,:)eXZ 
x>0 


S.t. 


Zj  =  zj,\/j e  Vt,i  =  l,...,r-l  (if r>l) 

<  Zj  g  {0,  l} ,  V/  g  Vr 

Zj£  [0,1], V/e  Vi,i  =  r  +  \,...,k  (if  r<k) 


where  the  values  zj  for  j e  f ,  /  =  1  - 1  in  stage  r>  1  are  retrieved  from  the  solution  to 

problems  SCIPA1,...,SCIPAr"1 ,  respectively.  Because  only  a  reduced  subset  of  (non- 
fixed)  0-1  variables  are  restricted  to  be  integer  at  each  stage,  r,  we  expect  each  of  the 
SCIPAr  models  to  solve  more  efficiently  than  original  SCIPA. 

In  particular,  our  implementation  begins  by  relaxing  the  binary  constraints  for  all 
the  variables  in  z  except  those  associated  with  period  y  =  1 .  SCIPA1  makes  it  possible  to 
easily  obtain  a  “what-to-do-fi rst”  solution.  These  binary  variables  are  then  fixed  at  the 
second  stage.  In  SCIPA2 ,  only  those  variables  associated  with  the  second  period 
( y  =  2)  are  deemed  integer.  This  cascade  process  is  followed  until  the  variables  for  the 

last  period,  y  =  \Y\ ,  are  set  to  integer  values. 

In  short,  our  model  (SCIPA)  is  divided  into  A:  =  |f|  sub-problems  that  need  to  be 
solved  in  sequence: 
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SCIPAy :  mm  f(a) 

(1)  to  (44) 

APROCayi  >  0 
ARETay  >  0 
SRETy  >  0 

AP  =  AP  ,\/ae  A,ie  I  ; y  =  1,...,  v-1  (if  v  >  1) 

ay  i  ay  i  a  y  y  \  s  / 

S.t.< 

SPROC  .  =  SPROC  ,  ,Vje  S,pe  P;VqeQ  v'  =  1,...,  v -1  (if  v  >  1) 
},Vas4is/, 

SPROC^e  {0,1}, Vss  S.pe  P,;Vq s 
-IP,  s|(l,l],  W(E  4is  /„;y  =y  +  l,...,|7|  (ify<|F|) 

SPROC  £  [0,1], Vs£  s,p*  e  =  v  + 1 . |r|  (if y  <  |r|) 

Apy  */  *py  11  1  ' 


If  F*  (SCIPA)  is  allowed  to  denote  the  optimal  objective  function  value  for  our 
original  model,  and  F(SCIPA)  and  F(SCIPA)  are  also  allowed  to  denote  a  lower 

bound  and  an  upper  bound  on  that  solution,  respectively,  the  F&R  algorithm  is  as 
follows: 

F&R  (SCIPA):  Fix-and-Relax  Algorithm  for  model  SCIPA 

Input:  Partition^, ...,Vk,  where  k  =  |7| ,  and  each  V y  contains  exactly  all 

the  binary  variables  associated  with  period  y: 

Vy=  {triplets  (a,  y,  i)  for  A Pavi  variables}  U  {four-uplas  (s,  p,  y,  q) 
for  SPROCSpyq  variables} 

Step  1:  Set  v=l  and  solve  SCIPA3' 

If  SCIPA1'  is  infeasible,  STOP:  “Problem  SCIPA  is  infeasible”. 
Otherwise,  set  F(SCIPA)  =  V* (SCIPA  ) . 

Step  2:  If y=k,  set  F(SCIPA)  =  F(SCIPAa)  and  STOP:  “Problem  SCIPA 
is  feasible”. 

Otherwise,  increase  y  by  1 . 
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Step  3:  Solve  SCIPA . 

If  SCIPA’  is  infeasible,  STOP:  “Problem  SCIPA  status  is 
unknown”. 

Otherwise,  go  back  to  Step  2 

Output:  SCIPA  status  (“Infeasible”,  “Feasible”  or  “Unknown”).  If  status 
is  “Feasible”,  U(SCIPA)  and  U(SCIPA)  are  a  lower  and  an 
upper  bound,  respectively,  on  the  optimal  solution  to  SCIPA. 

As  indicated  in  Step  3,  F&R(SCIPA)  has  the  potential  to  fail.  This  may  occur  if 
(SCIPA1)  is  feasible  but,  at  some  stagey  >1  ,  the  associated  problem  (SCIPAy)  becomes 
infeasible.  In  this  situation,  F&R(SCIPA)  is  unable  to  recognize  if  the  infeasibility  is  due 
to  the  fact  that:  (a)  SCIPA  is  actually  integer-infeasible  (but  continuous-feasible),  or  (b) 
(SCIPA)  is  integer-feasible,  but  the  cascade  fixing  procedure,  which  works  with 
estimates  of  the  true  optimal  values  of  the  variables,  makes  (SCIPAy)  infeasible. 

In  our  computational  experience,  the  later  problem  never  occurred,  but  it  if  did, 
alternative  versions  of  this  algorithm  may  be  implemented  that  overcome  this  difficulty 
(e.g.,  Escudero  and  Sahneron  [2002]). 

Notice  also  that  F&R(SCIPA)  yields  a  relative  gap  equal  to 

(U(SCIPA)  -F(SCIPA))/U  (SCIPA) .  Enhancements  of  the  algorithm  (e.g.,  Escudero 

and  Sahneron  [2002])  can  deal  with  the  situation  where  this  gap  is  too  big.  One  technique 
consists  of  stepping  back  and  grouping  multiple  stages  into  a  single  one.  Eventually,  if 
the  gap  discrepancies  continue,  the  F&R  method  becomes  a  single-stage  process  solving 
the  original  SCIPA. 
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III.  RESULTS 


This  chapter  analyzes  the  results  extracted  from  multiple  runs  of  CIPA  for 
different  test  cases. 

GAMS  [Brooke  et  al.  (1998)]  (version  2. 0.8. 3  with  Revision  117  module) 
incorporating  the  CPLEX  solver  [GAMS-CPLEX  (2002)]  (version  6.6.1)  solves  the 
MILP.  Computations  are  from  a  Dell  Computer  Precision  340  Pentium-4,  2  GHz  desktop 
computer  with  1  GB  of  random  access  memory. 

Two  main  data  sets  are  used  for  testing  purposes.  The  first  is  from  Baran  [2000], 
and  called  “Baseline  scenario  #1.”  The  second,  provided  by  N81,  is  called  “Baseline 
scenario  #2.” 

Baseline  scenario  #1  consists  of  30  aircraft  types,  12  different  air  missions,  13 
plants,  45  ship  classes,  and  17  different  ship  missions.  The  original  data  set  is  modified 
slightly  to  accommodate  set-aside  budget  and  labor  data. 

Baseline  scenario  #2  consists  of  38  aircrafts  types,  10  different  air  missions,  7 
plants,  46  ship  classes,  and  10  different  ship  missions.  This  data  set  has  not  been  used  in 
any  previous  thesis  research. 

We  create  24  different  cases  (1  0,. . .,1  3,. . .,6  0,. . .,6  3)  from  each  baseline 
scenario.  Because  we  use  the  same  notation  to  represent  the  same  type  of  excursion  from 
each  baseline  scenario,  results  are  presented  in  separate  tables  for  each  scenario.  In 
particular,  excursions  are  created  as  follows  (in  either  scenario):  Mission  requirements 
are  increased  10%,  25%  and  are  decreased  15%  to  create  cases  11,  1_2  and  13, 
respectively.  In  addition  to  the  mission  requirement  increment  (MRI),  the  budget  is 
decreased  20%  for  each  case,  yielding  cases  2_0  through  2_3.  An  individual  yearly 
budget  (IB)  option  is  added  to  cases  in  group  1,  creating  cases  3  0  through  3  3.  The 
same  modification  is  made  to  group  2  in  order  to  create  cases  4  0  through  4  3.  Groups  5 
(cases  5_0  through  5_3)  and  6  (cases  6_0  through  6_3)  are  created  by  removing  the 
cumulative  budget  (CB)  option  from  groups  3  and  4,  respectively. 
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Each  run  for  cases  under  scenario  #1  contains  about  14,000  continuous  variables 
and  4,300  discrete  variables.  Each  run  for  cases  under  scenario  #2  contains  about  15,000 
continuous  variables  and  5,100  discrete  variables. 

In  order  to  assess  the  efficiency  of  F&R,  results  are  compared  with  two  other 
methods:  The  first  method  is  a  branch-and-bound  (B&B)  [e.g.,  Wolsey,  1998],  as 
implemented  using  default  settings  in  GAMS-CPLEX  [2002].  In  order  to  be 
comprehensive,  a  post-rounding  process  of  the  solution  to  SCIPA  is  incorporated  in  order 
to  attain  feasibility  of  the  original  CIPA  model.  In  essence,  the  B&B  method  can  be 
viewed  as  a  F&R  instance  with  one  unique  stage,  i.e.,  a  single  partition  spanning  all  the 
integer  variables  of  the  problem.  The  second  method  is  a  customized  local-search 
heuristic  [Sahneron  et  ah,  2002]  that  typically  finds  acceptable  solutions  quickly. 

Tables  are  presented  in  a  compact  way  in  the  Appendix.  For  simplicity  of 
exposition,  in  this  chapter  we  divide  these  tables  by  the  following  measures  of 
effectiveness:  Upper  Bound  (UB),  Lower  Bound  (LB)  and  gap. 

A.  COMPARING  UPPER  BOUNDS 

Tables  1  and  2  show  the  results  for  the  two  scenarios  and  the  best  feasible 
solution  (i.e.,  the  best  UB)  obtained  using  the  three  methods. 
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Baseline  Scenario  #1 


CASE 

IB 

CB 

MRI 

Bl 

RB(Heur) 

(1) 

BB(B  &B) 
(2) 

UB  (F&R) 
(3) 

Time  for 
UB(B&B) 
(sec)  (2)(4) 

Time  for 
UB(F&R) 
(sec)  (3)(4) 

1_0 

* 

123536 

129203 

144367 

177.83 

434.90 

1  _1 

* 

10% 

730261 

753120 

751905 

254.06 

330.78 

1_2 

* 

25% 

1844213 

1917775 

1911800 

24.02 

367.48 

1_3 

* 

-15% 

75639 

98967 

126558 

717.84 

797.38 

2_0 

* 

-20% 

124127 

161590 

138476 

617.88 

565.15 

2_1 

UP 

10% 

-20% 

838422 

799877 

837392 

271 .89 

464.36 

2_2 

* 

25% 

-20% 

2240723 

2273176 

2305869 

74.48 

541.94 

2_3 

* 

-15% 

-20% 

76390 

100936 

124498 

438.51 

799.33 

3_0 

* 

* 

177610 

No 

181272 

1201.14 

1419.98 

3_1 

i 

* 

10% 

786317 

861380 

812775 

172.67 

473.68 

3_2 

* 

* 

25% 

1914442 

2001053 

2015089 

45.22 

283.13 

3_3 

* 

* 

-15% 

117154 

No 

115137 

1200.95 

1454.37 

4_0 

UP 

UP 

-20% 

144655 

200819 

185948 

1163.34 

684.35 

4_1 

UP 

* 

10% 

-20% 

880290 

859989 

929985 

717.51 

633.40 

4_2 

* 

* 

25% 

-20% 

2300018 

2289095 

2338256 

78.39 

492.09 

4_3 

UP 

UP 

-15% 

-20% 

95100 

No 

108379 

1200.98 

1038.1  1 

5_0 

UP 

174085 

No 

184313 

1200.95 

1332.08 

5_1 

UP 

10% 

784407 

832010 

833147 

139.00 

510.13 

5_2 

UP 

25% 

1915798 

1977670 

2005881 

36.01 

402.40 

5_3 

UP 

-15% 

120080 

No 

105776 

1200.95 

1621.90 

6_0 

UP 

-20% 

142908 

239287 

187633 

1205.84 

720.35 

6_1 

* 

10% 

-20% 

851852 

853657 

784293 

104.83 

824.42 

B_2 

* 

25% 

-20% 

2178379 

2028560 

2097786 

36.78 

389.13 

6_3 

UP 

-15% 

-20% 

9801  1 

No 

112660 

1200.91 

1031.47 

(1)  Heuristic  CPU  time  is  always  about  1  minute. 

(2)  B&B  runs  with  a  20  minute  time  limit.  "No"  indicates  the  case  that  B8B  cannot  obtain  a  feasible  solution  within  the 
allotted  time. 

(3)  Each  stage  of  F&R  obtains  the  solution  using  a  3  minute  time  limit  with  the  total  time  for  all  stages  shown. 

(4)  Includes  a  (negligible)  time  to  round  the  solution  for  SRET,  ARET,  APROC  variables. 

IB:  Individual  Budget;  CB:  Cumulative  Budget;  MRI:  Mission  Requirement  Increment;  Bl:  Budget  Increment 

Table  1.  Upper  Bound  (UB)  and  Computational  Run  Time  for  All  Methods  and  All 

Excursions  from  Baseline  Scenario  #1. 
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Baseline  Scenario  U2 


CASE 

IB 

CB 

MRI 

Bl 

UB(Heur) 

(1) 

UB(B&B) 

P) 

UB  (F&R) 
(3) 

Time  for 

(sec)  (2) (4) 

Time  for 
UB  (F&R) 
(sec)  (3) (4) 

1_0 

* 

1 526377 

No 

1538233 

2400.00 

2087.24 

1  _1 

* 

10% 

1 879666 

No 

1875393 

2400.00 

1551.54 

1_2 

25% 

2495731 

2511077 

2490478 

391.45 

1649.67 

1_3 

i 

-15% 

1072505 

No 

1068038 

2400.00 

1467.56 

2_0 

i 

-20% 

1 53 1 306 

1542857 

1538894 

502.43 

2220.30 

2_1 

i 

10% 

-20% 

1884332 

1880542 

1889328 

402.32 

1694.45 

2_2 

i 

25% 

-20% 

2508233 

No 

2499054 

2400.00 

1237.33 

2_3 

* 

-15% 

-20% 

1075544 

No 

1067096 

2400.00 

1741.88 

3_0 

* 

* 

1536413 

1528043 

1536010 

518.43 

1426.44 

3_1 

* 

* 

10% 

1889014 

1899922 

1878403 

387.60 

1109.18 

3_2 

* 

* 

25% 

2508594 

2536668 

2495177 

517.42 

1393.43 

3_3 

* 

* 

-15% 

1079994 

1075390 

1067605 

499.40 

1584.10 

4_0 

* 

* 

-20% 

1544200 

No 

1601  896 

2400.00 

1685.10 

4_1 

* 

* 

10% 

-20% 

1895842 

No 

1884112 

2400.00 

1316.77 

4_2 

* 

* 

25% 

-20% 

2520067 

No 

2555419 

2400.00 

1430.99 

4_3 

* 

* 

-15% 

-20% 

1086630 

1084394 

1069629 

492.29 

1494.43 

5_0 

* 

1529508 

No 

1523776 

2400.00 

1454.77 

5_1 

* 

10% 

1880137 

1883598 

1878234 

632.81 

1424.91 

5  _2 

* 

25% 

2501399 

2499434 

2497334 

383.98 

1400.31 

5_3 

* 

-15% 

1072848 

1074187 

1067744 

423.76 

1375.98 

6_0 

* 

-20% 

1544863 

No 

1536329 

2400.00 

1357.54 

6_1 

* 

10% 

-20% 

1896190 

1886009 

1889655 

491.33 

1478.48 

6_2 

* 

25% 

-20% 

2521073 

2520460 

2506490 

583.73 

1632.49 

6_3 

* 

-15% 

-20% 

1087255 

No 

1072530 

2400.00 

1376.93 

(1)  Heuristic  CPU  time  is  always  about  3  minutes. 

(2)  B&B  runs  with  a  40  minute  time  limit.  "No"  indicates  the  case  that  B&B  cannot  obtain  a  feasible  solution  within  the 
allotted  time. 

P)  Each  stage  of  F&R  obtains  the  solution  using  a  3  minute  time  limit  with  the  total  time  for  all  stages  shown. 

(4)  Includes  a  (negligible)  time  to  round  the  solution  for  SRET,  ARET,  APROC  variables. 

IB:  Individual  Budget;  CB:  Cumulative  Budget;  MRI:  Mission  Reguirement  Increment;  Bl:  Budget  Increment 

Table  2.  UB  and  Computational  Run  Time  for  All  Methods  and  All  Excursions  from 

Baseline  Scenario  #2. 


There  are  several  cases  that  cannot  be  solved  with  B&B,  or  even  obtain  a  feasible 
solution  within  20  minutes  (scenario  #1)  and  40  minutes  (scenario  #2).  Moreover,  in  most 
of  these  cases,  no  solution  at  all  is  obtained  even  if  this  method  is  executed  for  several 
hours.  By  employing  the  F&R  methodology,  feasible  solutions  are  obtained  within, 
approximately,  10-20  minutes  for  scenario  1,  and  20-40  minutes  for  scenario  #2. 

The  upper  bound  provided  by  the  B&B  (if  any)  and  F&R  methods  is  compared  to 
that  of  a  customized  heuristic  for  the  problem.  In  most  cases  under  scenario  #1,  the 
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heuristic  bound  is  superior.  For  scenario  #2,  the  F&R  bound  is,  in  general,  the  best, 
although  both  the  heuristic  and  B&B  bounds  are  very  close  to  that  value. 

It  is  verified  by  inspection  that  the  trade-off  between  computational  time  and 
solution  value  achieved  favors  the  use  of  the  heuristic  method  to  compute  such  solution. 
However,  it  cannot  be  ruled  out  that,  for  other  scenarios  or  future  specifications  of  the 
problem,  the  heuristic  will  not  behave  as  well  as  for  scenarios  #1  and  #2  which  served  as 
test-cases  for  its  development.  In  such  a  situation,  and  with  B&B  not  being 
computationally  affordable,  a  general  “quasi-exact”  approach,  such  as  F&R,  may  still  be 
needed. 

B.  COMPARING  LOWER  BOUNDS 

Solving  the  linear  programming  (LP)  relaxation  of  CIPA  (or  SCIPA)  takes  only 
about  one  minute  for  both  scenarios  #1  or  #2.  When  using  B&B,  a  better  LB  can  be 
obtained  by  inspecting  the  B&B  tree  and  selecting  the  best,  least  cost,  active  node.  For 
the  F&R  method,  the  best  LB  is  provided  by  the  optimal  solution  to  stage  1.  These  three 
bounds  are  illustrated  in  Tables  3  and  4. 

As  Tables  3  and  4  show,  the  three  LBs  look  very  similar  in  all  cases,  and  in  fact 
coincide  for  the  LP  relaxation  and  F&R. 
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Baseline  Scenario  #1 


CASE 

IB 

CB 

MRI 

Bl 

LP 

Relaxation 

LB(B&B) 

LB(F&R) 

1_0 

k 

110578 

119628 

119578 

1  _1 

k 

10% 

718213 

71  8359 

718213 

1_2 

* 

25% 

1774185 

1774201 

1774185 

1_3 

k 

-15% 

72854 

73000 

72954 

2_0 

k 

-20% 

113578 

119673 

119578 

2_1 

k 

10% 

-20% 

734973 

742579 

734973 

2_2 

k 

25% 

-20% 

2049108 

2068374 

2049108 

2_3 

k 

-15% 

-20% 

72954 

73000 

72954 

3_0 

* 

k 

150432 

150432 

150432 

3_1 

* 

k 

10% 

750281 

751046 

750281 

3_2 

k 

k 

25% 

1828030 

1828053 

1828030 

3_3 

k 

k 

-15% 

103747 

103747 

103747 

4_0 

k 

k 

-20% 

135443 

135566 

135443 

4_1 

* 

k 

10% 

-20% 

751577 

762405 

751577 

4_2 

* 

k 

25% 

-20% 

2070090 

2088395 

2070090 

4_3 

* 

k 

-15% 

-20% 

88766 

88766 

88766 

5_0 

k 

150432 

150432 

150432 

5_1 

k 

750281 

750637 

750281 

5_2 

k 

1828030 

1828056 

1828030 

5_3 

k 

-15% 

103747 

103747 

103747 

6_0 

k 

-20% 

135443 

135566 

135443 

6_1 

k 

-20% 

741844 

744259 

741844 

6_2 

k 

-20% 

1852129 

1855145 

1852129 

6_3 

k 

-15% 

-20% 

88766 

88766 

88766 

IB:  Individual  Budget;  CB:  Cumulative  Budget;  MRI:  Mission  Requirement  Increment;  Bl:  Budget  Increment 

Table  3.  Lower  Bound  (LB)  for  Baseline  Scenario  #1. 
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Baseline  Scenario  #2 


CASE 

IB 

CB 

MRI 

Bl 

LP 

Relaxation 

LB(B&B) 

LB(F&R) 

1_0 

* 

1520923 

1520923 

1520923 

1  _1 

* 

10% 

1874149 

1874149 

1874149 

1_2 

* 

25% 

2487142 

2487142 

2487142 

1_3 

* 

-15% 

1066127 

1066127 

1066127 

■ 

H 

-20% 

1520923 

1521029 

1520923 

2_1 

* 

10% 

-20% 

1874149 

1874149 

1874149 

2_2 

* 

25% 

-20% 

2487142 

2487142 

2487142 

2_3 

* 

-15% 

-20% 

1066127 

1066127 

1066127 

H 

H 

1520923 

1520923 

1520923 

1874149 

1874149 

1874149 

2487142 

2487193 

2487142 

3_3 

* 

* 

-15% 

1066127 

1066127 

1066127 

ESI 

H 

H 

-20% 

1520923 

1520923 

1520923 

4_1 

* 

* 

10% 

-20% 

1874149 

1874149 

1874149 

4_2 

* 

* 

25% 

-20% 

2487142 

2487142 

2487142 

4_3 

* 

* 

-15% 

-20% 

1066127 

1066127 

1066127 

5_0 

* 

1520923 

1520923 

1520923 

5_1 

* 

1874149 

1874149 

1874149 

5_2 

* 

2487142 

2487143 

2487142 

5_3 

* 

-15% 

1066127 

1066127 

1066127 

m 

H 

-20% 

1520923 

1520923 

1520923 

6_1 

0 

-20% 

1874149 

1874149 

1874149 

H 

-20% 

2487142 

2487142 

2487142 

6_3 

* 

-15% 

-20% 

1066127 

1066127 

1066127 

IB:  Individual  Budget;  CB:  Cumulative  Budget;  MRI:  Mission  Requirement  Increment;  Bl:  Budget  Increment 

Table  4.  LB  for  Baseline  Scenario  #2. 

C.  GAP  COMPARISON 

In  order  to  assess  the  quality  of  the  solution  obtained,  we  define  a  relative  gap  for 
jjb _ m 

each  method  as  the  ratio  - ,  where  UB  and  LB  refer  to  the  upper  and  lower 

LB 

bounds,  respectively,  provided  by  the  method.  For  the  heuristic  method,  we  take  the  LP 
relaxation  as  a  lower  bound. 

Overall,  the  heuristic  methodology  yields  the  best  solution  among  the  three 
methods  for  scenario  #1  (see  Table  5)  and  provides  the  solution  within  acceptable 
tolerance  for  scenario  #2  (see  Table  6).  This  is  not  surprising  since  the  heuristic  solver 
was  developed  and  tested  using  these  scenarios  as  training  cases. 

As  shown  on  Table  6  for  base  line  scenario  #2,  F&R  is  more  accurate  than  the 

other  methods.  In  particular,  B&B  cannot  reach  a  feasible  solution  within  a  given  time  in 
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many  instances  of  this  scenario.  Although  F&R  is  the  overall  winner  for  scenario  #2,  the 
heuristic  solution  is  as  almost  the  same  quality  as  the  F&R  solution. 


Baseline  Scenario  #1 


CASE 

IB 

CB 

MRI 

Bl 

Gap 

(Heur) 

Gap 

(B&B) 

(i) 

Gap 

(F&R) 

1_0 

* 

3.31% 

8.00% 

20.73% 

1_1 

* 

10% 

1.68% 

4 .84% 

4.69% 

1_2 

* 

25% 

3.95% 

8.09% 

7.76% 

1_3 

* 

-15% 

3.68% 

35.57% 

73.48% 

2_0 

* 

-20% 

3.80% 

35.03% 

15.80% 

2_1 

* 

10% 

-20% 

14.08% 

7  72% 

13.94% 

2  2 

* 

25% 

-20% 

9.35% 

9.90% 

12.53% 

2_3 

* 

-15% 

-20% 

4.71% 

38.27% 

70.65% 

3_0 

* 

* 

18.07% 

? 

20.50% 

3_1 

* 

* 

10% 

4.80% 

14.69% 

8.33% 

3_2 

* 

* 

25% 

4.73% 

9.46% 

10.23% 

3_3 

* 

* 

-15% 

12.92% 

? 

10.98% 

4_0 

* 

* 

-20% 

6.80% 

48.13% 

37.29% 

4_1 

* 

* 

10% 

-20% 

17.13% 

12.80% 

23.74% 

4  2 

* 

* 

25% 

-20% 

11.11% 

9.61% 

12.95% 

4_3 

* 

* 

-15% 

-20% 

7.14% 

? 

22.10% 

5_0 

* 

15.72% 

? 

22.52% 

5  1 

* 

10% 

4.55% 

10.84% 

1 1 .04% 

5_2 

* 

25% 

4.80% 

8.18% 

9.73% 

5_3 

* 

-15% 

15.74% 

? 

1 .96% 

6_0 

* 

-20% 

5.51% 

76.51% 

38.53% 

6_1 

* 

10% 

-20% 

14.83% 

14.70% 

5.72% 

6  2 

* 

25% 

-20% 

17.61  % 

9.35% 

13.26% 

6_3 

* 

-15% 

-20% 

10.42% 

? 

26.92% 

(1)  The  MlPto  SCIPA  solution  is  within  10%  tolerance.  However  after  rounding  for  the  SRETT,  ARET,  APROC  variables  the 
gap  may  increase  to  more  than  10%. 

IB:  Individual  Budget;  CB:  Cumulative  Budget;  MRI:  Mission  Reguirement  Increment;  Bl:  Budget  Increment 

Table  5.  Gap  between  LB  and  UB  for  Baseline  Scenario  #1. 

The  unexpected  results  in  Tables  5  and  6  are  the  cases  in  which  the  gap  is  bigger 
than  10%  for  F&R.  We  next  discuss  the  reason  for  these  large  gaps  by  observing  the 
objective  value  of  case  13  in  base  line  scenario  #1  throughout  the  evolution  of  the  F&R 
algorithm. 
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Baseline  Scenario  #2 


CASE 

IB 

CB 

MRI 

Bl 

Gap 

(Heur) 

Gap 

(B&B) 

(1) 

Gap 

(F&R) 

1_0 

* 

0.36% 

? 

1.14% 

1_1 

* 

10% 

0.29% 

? 

0.07% 

1_2 

* 

25% 

0.35% 

0.96% 

0.13% 

1_3 

* 

-15% 

0.60% 

? 

0.18% 

2_0 

* 

-20% 

0.63% 

1 .44% 

1.18% 

2_1 

* 

10% 

-20% 

0.54% 

0.34% 

0.81% 

2_2 

* 

25% 

-20% 

0.85% 

? 

0.48% 

2_3 

* 

-15% 

-20% 

0.88% 

? 

0.09% 

3_0 

* 

1.02% 

0.47% 

0.99% 

3_1 

* 

* 

10% 

0.79% 

1 .38% 

0.23% 

3_2 

* 

25% 

0.86% 

1 .99% 

0.32% 

3_3 

* 

* 

-15% 

1.30% 

0.87% 

0.14% 

4_0 

* 

-20% 

1.53% 

? 

5.32% 

4_1 

* 

* 

10% 

-20% 

1.16% 

? 

0.53% 

4_2 

* 

25% 

-20% 

1.32% 

? 

2.75% 

4_3 

* 

* 

-15% 

-20% 

1.92% 

1.71% 

0.33% 

5_0 

T 

0.56% 

? 

0.19% 

5_1 

10% 

0.32% 

0.50% 

0.22% 

5_2 

* 

25% 

0.57% 

0.49% 

0.41% 

5_3 

* 

-15% 

0.63% 

0.76% 

0.15% 

6_0 

-20% 

1.57% 

? 

1.01% 

6_1 

10% 

-20% 

1 .1 8% 

0.63% 

0.83% 

6_2 

* 

25% 

-20% 

1.36% 

1 .34% 

0.78% 

6_3 

-15% 

-20% 

1.98% 

? 

0.60% 

(1)  The  MlPto  SCIPA  solution  is  within  10%  tolerance.  However  after  rounding  for  the  SRETT,  ARET,  APROC  variables  the 
gap  may  increase  to  more  than  10%. 

IB:  Individual  Budget;  CB:  Cumulative  Budget;  MRI:  Mission  Requirement  Increment;  Bl:  Budget  Increment 

Table  6.  Gap  between  LB  and  UB  for  Baseline  Scenario  #2. 


As  Figure  1  and  Table  7  show,  the  gap  for  the  F&R  solutions  is  within  acceptable 
limits.  The  gap  increases  suddenly  in  the  post-rounding  process,  which  is  used  to  attain 
integer  feasibility  for  the  original  CIPA  model.  So,  the  F&R  procedure  we  have 
developed  is  actually  near-optimal  if  we  consider  the  model  that  it  has  been  applied  to: 
SCIPA.  The  gap  for  the  final  CIPA  can  attributed  to  the  rounding  process. 
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Figure  1.  Stage-by-Stage  F&R  Gap  for  Case  1_3  in  Scenario  #1. 

Case  1  3  in  Baseline  Scenario  #1 


F&R 

Stage 

Obj  .Value 

Gap 

Time(sec) 

FY06 

72954 

0.00% 

7.72 

FY07 

72954 

0.00% 

30.06 

FY08 

72954 

0.00% 

23.91 

FY09 

73015 

0.08% 

41.89 

FYTO 

73066 

0.15% 

20.69 

FY1 1 

73107 

0.21% 

43.53 

FYT2 

73107 

0.21% 

55.17 

FY13 

73125 

0.23% 

21.22 

FY14 

73214 

0.36% 

47.81 

FY'15 

73214 

0.36% 

58.19 

FY16 

73214 

0.36% 

41.47 

FY17 

73265 

0.43% 

49.62 

FY18 

73539 

0.80% 

23.59 

F  Y 1 9 

73590 

0.87% 

37.39 

FY20 

73590 

0.87% 

33.81 

FY21 

73684 

1.00% 

26.31 

FY22 

73920 

1.33% 

52.59 

FY23 

73920 

1.33% 

17.33 

FY24 

74639 

2.31% 

39.33 

FY25 

76464 

4.81% 

17.78 

FY26 

76475 

4.83% 

34.94 

FY27 

76929 

5.45% 

33.27 

FY28 

77052 

5.62% 

15.17 

FY29 

77155 

5.76% 

9.33 

FY30 

77374 

6.06% 

12.06 

Rounding 

126558 

73.48%i 

3.20 

Total  797.38 


B&B 

Stage 

Obj.  Value 

Gap 

B&B 

75618 

3.59% 

Rounding 

98967 

35.57% 

Table  7.  Stage-by-Stage  F&R  Gap  and  B&B  Gap  for  Case  13  in  Scenario  #1. 

36 


Table  7  raises  one  more  interesting  issue.  When  F&R  solves  one  stage,  it  can 
yield  an  integer  solution  for  some  subsequent  stages.  This  suggests  that  part  of  the  time  to 
solve  the  problem  for  these  stages  can  be  saved.  For  instance,  the  objective  function  value 
for  FY06  does  not  change  for  the  next  two  steps.  If,  after  solving  FY06,  the  solution  for 
variables  indexed  by  FY07  and  FY08  is  verified  to  already  be  integer,  the  solver  can 
continue  to  find  a  feasible  solution  by  skipping  these  stages.  This  is  also  valid  for  FY1 1- 
12,  FY14-16,  FY19-20,  and  FY22-23.  If  the  strategy  for  the  case  shown  in  Table  7  had 
been  implemented,  the  total  run  time  would  become  537.44  seconds  instead  of  797.38 
seconds. 
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IV.  END-EFFECTS  ANALYSIS 


Working  with  a  finite  horizon  is  a  simplification  driven  by  problem  complexity 
and  limited  knowledge  of  data.  Unfortunately,  in  many  cases,  using  an  artificial  finite 
horizon  adversely  influences  the  optimal  decisions,  referred  to  as  end-effects  [e.g., 
Walker  et  ah,  1995],  CIPA  optimizes  the  problem  over  a  30-year  planning  horizon,  but 
the  actual  procurement  and  retirement  of  U.S.  Navy  assets  should  extend  beyond  this 
point.  End-effects  that  emerge  because  of  the  mentioned  situation  arise  especially  due  to 
two  reasons:  (1)  because  no  future  missions  are  visualized  after  the  last  year,  and  (2) 
because  the  cost  and  labor  structure  of  some  platfonns  hinder  spending  money  or  labor 
for  deliveries. 

To  overcome  this  problem,  the  concepts  of  “set  aside  budget”  (for  ships  and 
aircraft)  and  “set  aside  labor”  (for  ships)  are  added  into  the  structure  of  CIPA.  In  this 
context,  the  planner  may  specify  maximum  amounts  of  these  categories  to  be  set  aside  for 
(undecided)  procurements  in  years  beyond  the  plan’s  scope.  The  maximum  labor  to  be  set 
aside  is  specified  by  plant  and  year  [Sahneron  et  ah,  2002].  Also,  a  consistent  relationship 
between  the  set-aside  budget  and  set-aside  labor  is  enforced. 

The  change  has  been  implemented  in  the  CIPA  model  presented  in  Chapter  III 
(Equations  (18),  (23)  and  (26)).  In  this  chapter,  how  the  model  behaves  with  and  without 
end  effects  is  explored. 

A.  OVERALL  MISSION  EFFECTIVENESS  WITH  AND  WITHOUT  SET- 

ASIDES 

CIPA  aims  to  minimize  budget,  industry,  and  mission  requirement  violations.  One 
way  to  reduce  minimum  budget  penalties  is  to  spend  money  procuring  platforms  without 
mission  requirements. 
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AAW  effectiveness  for  aircrafts 


Figure  2.  Anti  Aircraft  Warfare  (AAW)  Effectiveness  for  Aircrafts  Incorporating 
Mission  Requirements,  the  Case  with  and  without  Set-Aside  Concepts. 


For  example,  Figure  2  illustrates  that  the  model  without  set-aside  concepts 
minimizes  the  penalties  by  spending  more  money  to  procure  and  maintain  aircraft  assets 
even  if  they  are  not  needed  to  satisfy  the  given  mission  requirements.  In  contrast,  the 
model  with  set-aside  concepts  follows  plausible  near  mission  requirements.  This  situation 
is  likely  if,  during  some  years  of  the  planning  horizon,  there  is  a  budget  surplus  mismatch 
with  respect  to  mission  or  platform  requirements.  Allowing  set  asides  to  be  part  of  the 
budget  in  the  last  years  of  the  horizon  reduces  unnecessary  expenditures  in  those  years 
and  also  in  previous  years,  without  violating  yearly  and  cumulative  budget  limits.  (There 
are  other  missions  where  the  improvement  by  using  set-aside  is  not  as  noticeable  as  for 
AAW.) 

The  impact  of  end-effects  is  diminished  by  using  set  aside  budget  and  set  aside 
labor.  Without  this  feature,  the  results  may  be  misled  by  end-effects.  Consequently,  we 
are  getting  closer  to  reality  by  adding  the  set  aside  budget  and  set  aside  labor  concepts  to 
CIPA. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR 

FURTHER  STUDIES 

This  thesis  shows  that  F&R  is  a  reliable  methodology  for  solving  CIPA  MILPs. 
The  structure  of  this  MILP  lends  itself  to  a  time -based  decomposition,  which  is  the 
approach  used  for  implementing  F&R.  However,  F&R  offers  a  variety  of  alternatives  that 
may  work  even  better  than  the  one-period-per-stage  strategy.  This  thesis  opens  a  new 
door  for  further  studies  based  on  our  current  implementation  outcomes.  For  example, 
F&R  could  be  improved  by  grouping  variables  by  ship  class,  aircraft  type  or  multiple 
time  periods.  We  may  also  define  a  “worth”  for  each  integer  decision  variable  in  the 
problem,  associated  with  the  cost  of  the  platform,  or  its  overall  effectiveness,  or  both. 
Variables  would  be  grouped  according  to  their  worth  (with  the  most  valuable  variables 
grouped  in  early  stages),  and  the  same  F&R  methodology  can  be  applied. 

Another  enhancement  to  our  current  implementation  is  checking  variables’  status 
before  solving  further  stages.  This  would  allow  us  to  skip  consecutive  stages  that  have  all 
their  variables  already  set  to  integer  values,  which  in  turn  decreases  the  total  run  time  of 
F&R. 

Even  though  F&R  may  fail  to  converge  to  a  solution  within  10%  tolerance,  we 
demonstrate  that  in  our  test  cases  this  gap  can  be  attributed  to  the  rounding  process  after 
solving  the  SCIPA  model,  rather  than  the  F&R  technique  itself.  This  suggests  revising 
the  rounding  process,  or  implement  F&R  directly  for  CIPA,  instead  of  for  SCIPA. 

The  existing  heuristic  method  gives  good  results  for  most  of  the  cases  tested  in 
this  thesis.  However,  it  should  be  taken  into  consideration  that  these  cases  were  used  as 
training  cases  to  develop  the  heuristic  solver.  Thus,  we  cannot  guarantee  that  this  method 
can  always  be  successfully  applied  to  CIPA. 

We  analyze  the  so-called  “set  aside  budget”  and  “set  aside  labor”  concepts  for 
diminishing  end-effects.  The  incorporation  of  these  concepts  into  our  MILP  yields  more 
realistic  results  by  providing  effectiveness  levels  closer  to  mission  requirement  levels 
and,  in  general,  by  not  over-expending  during  the  planning  horizon.  Considering  the 


41 


possibility  of  set-aside  budget  and  set-aside  labor  allows  us  to  adjust  the  procurement 
plan  in  order  to  achieve  the  same  level  of  mission  accomplishment  avoiding  while 
unnecessary  over-expenditures. 

As  a  future  extension,  we  realize  that  maximum  amounts  of  the  set-aside  budget 
and  labor  might  be  part  of  the  optimization  decisions,  rather  than  input  data.  Further 
research  may  also  implement  infinite-horizon  linear  programming.  More  realistic  results 
may  be  obtained  by  using  either  approach. 

As  a  result  of  this  research,  we  recommend  using  the  existing  heuristic  solver 
first,  along  with  any  exact  lower  bound.  If  this  does  not  provide  a  solution  within  the 
desired  tolerance,  then  F&R  should  be  used  as  a  second  solver.  If  F&R  cannot  find  a 
solution  within  tolerance,  then  B&B  can  be  used. 

It  is  important  to  keep  focusing  on  the  modeling  aspects  of  the  problem  to  ensure 
that  it  meets  the  U.S.  Navy  needs  as  closely  as  possible,  as  well  as  developing  efficient 
solving  techniques.  We  think  that  techniques  using  decomposition  methods  and  taking 
advantage  of  the  problem  structure  are  the  most  promising  ones. 
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Baseline  Scenario  #1 


APPENDIX.  COMPARISON  OF  BASE  LINE  SCENARIOS 

INCORPORATING  WITH  HEURISTIC,  BRANCH-AND-BOUND 
(B&B)  AND  FIX-AND-RELAX  (F&R)  SOLVERS 
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Comparison  of  All  Methods  and  All  Excursions  from  Baseline  Scenario  #1. 
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Table  8 


(3)  Each  stage  of  F&R  obtains  the  solution  using  a  3  minute  time  limit  with  the  total  time  for  all  stages  shown. 

(4)  Includes  a  (negligible)  time  to  round  the  solution  for  SRET,  ARET,  APROC  variables. 

(5)  The  MIP  to  SCIPA  solution  is  within  10%  tolerance.  However  after  rounding  for  the  SRET,  ARET,  APROC  variables  the  gap  may  increase  to  more  than  1 0%. 
B:  Baseline  Scenario  #1;  IB:  Individual  Budget;  CB:  Cumulative  Budget;  MRI:  Mission  Requirement  Increment;  Bl:  Budget  Increment 


Baseline  Scenario  #2 


Table  9.  Comparison  of  All  Methods  and  All  Excursions  from  Baseline  Scenario  #2. 
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(3)  Each  stage  of  F&R  obtains  the  solution  using  a  3  minute  time  limit  with  total  timeforall  stages  shown. 

(4)  Includes  a  (negligible)  time  to  round  the  solution  for  SRET,  ARET,  APROC  variables. 

(5)  The  MIP  to  SCIPA  solution  is  within  10%  tolerance.  However  after  rounding  for  the  SRET,  ARET,  APROC  variables  the  gap  may  increase  to  more  than  1 0% 
B:  Baseline  Scenario  #2;  IB:  Individual  Budget;  CB:  Cumulative  Budget;  MRI  Mission  Requirement  Increment;  Bl:  Budget  Increment 
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